DOTL:10.1193 1/guihaia.gxzw201806002 


Phylogenetic and morphological profile of Cladophora fracta 
(Cladophorophyceae, Chlorophyta) from Karst spings, in 
northern China 


HU Bianfang ', JI Li’, CHEN Le *, FENG Jia *, SHI Shengli* 


(1. Department of Biology, Jinzhong University, Jinzhong 030600, Shanxi, China; 2.College of Environment and 
Safety, Taiyuan University of Science and Technology, Taiyuan 030024, China; 3. School of Pharmaceutical 
Science, Shanxi Medical University, Taiyuan 030001, China; 4. School of Life Science, Shanxi University, Taiyuan 
030006, China) 


Abstract: Cladophora fracta, a filamentous green macroalgal epiphyte on rhodoliths, is described 
from five karst springs in North China. Although Cladophora species frequently appear in karst 
system, their genetic diversity, biogeographical affinities and physiological properties have not 
been well investigated in these environments. The specific objectives of this study were to: 1) 
describe the habitat of the Cladophora-like algae form the five karst springs; 2) identify the thallus 
to species level based on a combination of morphological characteristics and molecular sequence; 
and 3) explore the morphological influence of habitat. To elucidate the biogeographical patterns in 
Cladophora, both morphological and molecular evidence were compared of Cladophora 
specimens across five study sites. Analyses of partial small subunit (SSU) and large subunit (LSU) 
genes revealed that the studied 50 Cladophora specimens were genetically identical species and a 
total of 13 ribotypes were detected. The molecular sequencing results indicated that the examined 
species was highly homologous with Cladophora vagabunda, though they shared few 
morphological features. The genus didn't form a monophyletic clade but in three different clades 
both in SSU and LSU trees. The microscopic structure was more consistent with that of C. fracta. 
The Cladophora from the five karst springs did not show significant variation in cell dimensions. 


However, the species exhibited larger cell diameters than those reported from lakes. In addition, 
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the rhizoid-like branches are only observed in two locations (XA and ST). Considering the 
morphological characteristics, we therefore hold our species as Cladophora fracta. 
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Cladophora is a large and common green macroalgal genus belonging to the Cladophoraceae 
in the Cladophorales, whose members consist of branched or unbranched uniseriate filaments, and 
absence of akinete. Species of the green algal genus Cladophora (Cladophorales, Chlorophyta) are 
widely distributed from marine to freshwater habitats worldwide (Leliaert et al, 2009) and extend 
into cold temperature and polar waters (Boedeker et al, 2012). Most species grow attached to 
rocky substrate by rhizoidal cells, but they also can form extensive free-floating masses in 


eutrophic waters. 


The taxonomy of Cladophora species has been problematic because of the wide degree of 
morphological variation in response of different environmental conditons (Hayakawa et al., 2012). 
The Cladophorales is polyphyletic with representatives in three main clades (Siphonocladus clade, 
Cladophora clade and Aegagropila clade) (Hanyuda et al, 2002; Leliaert et al, 2003; Yoshii et al, 
2004), and the fourth lineage (Okellyaceae including only one species) has been proposed as a 
sister to the three main clades (Leliaert et al, 2009; Boedeker et al, 2012). More recently, 
molecular phylogenetic studies based on SSU and LSU rDNA sequences have provided the insight 
into the relationships within the Cladophorales (Hanyuda et al, 2002; Leliaert et al, 2003, 2007; 
Boedeker et al, 2012; Ichihara et al, 2013). However, the systematic relationship of genera remain 
poorly understood. Molecular phylogeny revealed that the traditional family and genus level 
classifications did not reflect the phylogenetic relationships (Leliaert et al, 2007, 2009). The 
Cladophora species are distributed in all three lineages rather than form a monophyletic group. 

Genus Cladophora contains a heterogeneous group of species that are very hard to tell apart 
and classify, mainly because of the high morphological plasticity and cryptic diversity (Gestinari 
et al, 2010). Cladophora is one of the most species-rich genera of green macroalgae and the 
species of the genus are morphologically highly variable (van den Hoek & Chihara, 2000). It is 
also difficult to define stable taxonomic characteristics, since these characters are influenced by 
habitat, age and environmental conditions. However, studies on the habitat and ecology of genus 
Cladophora remain rare, especially about those in karst springs. 

In this study, we collected Cladophora-like filamentous green algae from five karst springs in 
northern China. We performed molecular phylogenetic analyses based on nuclear-encoded small 
subunit (SSU), partial large subunit (LSU) sequences and the combined SSU and LSU sequences, 
as well as morphological observation to reveal the diversity of these filamentous green algae from 
karst springs. The specific objectives of this study were to: 1) describe the habitat of the 
Cladophora-like algae form the five karst springs; 2) identify the thallus to species level based on 
a combination of morphological characteristics and molecular sequence; and 3) explore the 


morphological influence of habitat. 


1 Materials and methods 


1.1 Plant materials 


A field survey was conducted at five typical karst springs in Shanxi Province, China in June 
and July 2015 to identify the habitat of Cladophora. Figure 1 showed the sampling locations of 
Cladophora in Shanxi province. Water temperature and pH were measured on site with a 
pH/EC/TDS/Temperature Tester (HI98129, HHANNA Instruments Inc., Italia). Specimens of 
Cladophora were collected from the five streams and Table 1 gave detail information on the 
sampling sites. Voucher species were deposited at Herbarium of Jinzhong University (JZU) and 
the voucher information was also given in Table 1. Specimens used for the morphological studies 
were preserved in freshwater containing 4% formalin or kept alive in freshwater, and those used 
for the molecular studies were frozen at —-20 °C. External features of thalli were observed under 


an Olympus BX51 bright field microscope (Olympus Co., Tokyo, Japan). 


Table 1 Collection information and GenBank accession numbers of SSU and LSU sequences for Cladophora 


fracta populations 


Code Collection information Vocher Accession No. No. of 
SSU LSU haplotypes 
XA Xin’an Spring, Lucheng City, Shanxi JZU150702 KX421221 KX421233 2 
Province, North China (36.27 °N, JZU150706  KX421222 KX421234 


113.26 °E).Water temp. 15-16 °C, pH 
7.15-7.7. TDs 463 mg * L". Coll. Hu B F. 6 


July. 2015 
SQW Shuiquanwan Spring, Xinzhou City, Shanxi JZU150601 KX421215 KX421228 4 
Province, North China (38.51 °N, JZU150603  KX421216 KX421229 


113.29 °E).Water temp. 10.4-15.6 °C, pH JZU150604 KX421217 KX421230 
6.3-7.6. TDs 395 mg * L'. Coll. Shi S L.21 JZU150608 KX421218 - 


June. 2015 
ST Shen’tou Spring, Shuozhou City, Shanxi JZU150716  KX421219 KX421231 2 
Province, North China (39.39 °N, JZU150719  KX421220 KX421232 


112.61 °E). Water temp.14.6-17 °C, pH 
7.3. TDs 485 mg * L". Coll. Shi S L.3 July. 


2015 
JC Source of Nanlao Spring, Jinci Temple, JZU150722 KX421210 KX421223 3 
Taiyuan City, Shanxi Province, North JZU150723 KX421211 KX421224 


China(37.42 °N, 112.26 °E). Water temp. JZU150727  KX421212 KX421225 
11.6-16.2 °C, pH 6.5-7.3. TDs 552 
mg * L'. Coll. Hu B F.25 July. 2015 
LC Longci Spring, Linfen City, Shanxi JZU150736  KX421213 KX421226 2 
Province, North China (36.08 °N, JZU150740 KX421214 KX421227 
111.38 °E). Water temp. 15.1—18.4 °C, pH 
6.5-7.2. TDs 506 mg * L”. Coll. Shi S L.5 
July. 2015 
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Note: Location abbreviations were as given in Table 1. 
Fig. 1 Map of Shanxi Province showing the sampling locations of Cladophora fracta 

1.2 DNA methods 

Specimens were frozen in liquid nitrogen until use and the total DNA was extracted with 
Aqua-SPIN Plant gDNA Isolation Mini Kit (Watson Biotechnologies, Inc) following the 
manufacturer’s instructions. Molecular phylogenetic analyses were carried out based on 
nuclear-encoded small subunit (SSU) and partial large subunit (LSU), and both genes were 
combined in partitioned alignment. The primer pairs used for amplifying SSU were based upon 
Teng (2011) (18S rDNAF: 5’-AAT GGC TCG GTA AAT CAG TT-3’ and 18S rDNAR: 5’-AGT 
TGA TGA CTC GCG CTT AC-3’). LSU primers were based upon Leliaert et al (2007) (C1: 


5’-ACC CGC TGA ATT TAA GCA TATC-3’ and D2: 5’-TCC GTG TTT CAA GAC GG-3’). The 
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standard polymerase chain reaction (PCR) analyses were carried out on a thermocycler (MyCycler 
Thermal, BIO-RAD, USA). For SSU sequences, PCR amplification was carried out as follows: 
initial denaturation at 95 °C for 5 min, 35 cycles at95 °C for 30 s, 56.9 °C for 30S and 72 °C 
for 1 min 30 s, and a final extension at 72 °C for 7 min. For LSU, PCR amplification was carried 
out as follows: initial denaturation at 95 °C for 5 min, 35 cycles at 95 °C for 30s, 63.8 °C for 
30 s and 72 °C for 30 s, and a final extension at 72°C for 7 min. PCR product was purified with 
Gel Extraction Mini Kit (Watson Biotechnologies, Inc) according to the manufacturer’s 
recommendation for direct sequencing. The PCR products were commercially sequenced by 
Shanghai Personal Biotechnology Co, Ltd. Additonal taxa from which obtained from Genbank 
were shown on Figure 2. 
1.3 Sequence analyses 

The Clustal-X 2.0 software was used to align the sequences (Thompson et al, 1997), and the 
program jModeltest (Posada, 2008; Guindon & Gascuel, 2003) was used to determine parameters 
for the maximum likelihood and Bayesian analyses. For the SSU, the model was as follows: 
GTR+I+G distance model, portion of invariable sites = 0.592 0; gamma distribution = 0.486 0; 
base frequencies A = 0.246 7, C = 0.221 0, G = 0.286 9, T = 0.245 5; and rate matrix A-C = 0.783 
2, A-G = 2.434 9, A-T = 2.096 1, C-G = 0.613 3, C-T = 6.101 0. For the LSU gene, the model 
was as follows: GTR + G distance model, gamma distribution = 0.378 0; base frequencies A = 
0.208 6, C = 0.284 6, G = 0.305 6, T = 0.201 2; and rate matrix A-C = 0.842 7, A-G = 2.347 0, 
A-T = 1.569 8, C-G = 0.478 6, C-T = 4.404 0. For the combination of the genes, the model was 
as follows: for the SSU partition portion, the model was as follows: GTR+I+G distance model, 
portion of invariable sites = 0.592 0; gamma distribution = 0.486 0; base frequencies A = 0.246 7, 
C = 0.221 0, G = 0.286 9, T = 0.245 5; and rate matrix A-C = 0.783 2, A-G = 2.434 9, A-T = 
2.096 1, C-G = 0.613 3, C-T = 6.101 0. For the LSU partition, the model was as follows: GTR + 
G distance model, gamma distribution = 0.394 0; base frequencies A = 0.211 5, C = 0.285 0, G= 
0.299 5, T = 0.203 9; and rate matrix A-C = 0.900 2, A-G = 2.917 3, A-T = 1.856 0, C-G = 0.569 
4, C-T = 5.437 1. The Maximum Likelihood (ML) analyses were conducted using PhyML 3.0 
(Guindon & Gascuel 2003). The robustness of trees obtained from ML analyses was estimated 
using bootstrap resampling with 1000 replicates (Felsenstein, 1985). The Bayesian analyses of the 


combined data were conducted using MrBayes version 3.1.2 (Huelsenbeck et al, 1996; Ronquist & 
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Huelsenbeck, 2003). A MCMCMC (Metropolis-coupled Markov chain Monte Carlo) algorithm 
running four Markov chains simultaneously was used to estimate the posterior probability of the 
phylogenetic trees. The Markov chains were started from a random tree and run for 10 000 000 
generations sampling every 1 000 generations for a total of 10 000 samples each run. The first 2 
500 samples from each run were discarded as burn-in. The consensus tree was reconstructed after 


burn-in of 25% generations. 


2 Results 


2.1 Habitat 

Karst spring, part of a karst system, is ususlly the end of a cave system at the place where a 
river cave reaches the Earth’s surface. Shanxi Province is the most abundant and typical region of 
karst system in northern China, where the exposed karst region was 2.6x10* km’, occupying 
nearly 17.5% of the total area of the whole province (Han et al, 1993). There were 18 large karst 
springs with the average discharge being more than 1.0 m/s, including the five investigated karst 
springs. Water temperature of the five streams ranged from 10.4 to 18.6 °C in July (except that for 
SQW was in June). As expected the water temperature showed some seasonal and geographical 
fluctuations. The pH of the five collecting sites was more consistent over the sampling period, 
ranging from 6.3 to 7.7. Total dissolved solids (TDS) showed not much fluctuation over the 
sampling sites, varying from a minimum of 395 mg/L at SQW where the water temperature was 
the lowest to a maximum of 552 mg/L at JC where the water temperature was the highest. 
2.2 Morphological observations 

The morphological analyses were based on three specimens of each location. Plants form 
bright green, clustered, floating or mostly attached to the rock substrates, and measured at about 
10cm in height. Thalli are composed of irregularly branched, uniseriate filaments.,. Cells are 
almost pear shaped with irregular swelling. The main axes of the cells were about 250.00-466.70 
um in length, and the length/width ratio was highly variable, in the range of 2.63-6.99 in regions 
with few cell divisions. The apical cells were about 378.64-966.67 um in length, and the 
length/width ratio was also highly variable, in the range of 3.84-13.33. The rhizoid-like branches 
were only observed in XA and ST. Although there were morphological variances in the samples 


from five sites, when we carefully checked the diagnostic characters and compared to the previous 


reports (Hu & Wei, 2006), the other features were common in all individuals detected in this study. 
Therefore, the samples were identified as Cladophora fracta. 


Table 2 Morphological comparison of cells from samples with Cladophora fracta 


Code Main axes cells (n=30) Apical cells(n=30) Rhizoid-like 


branches(n=30) 


Length(um) L/D ratio Length(um) L/D ratio Diameter (um) 


XA 250.00-402.60 3.32-4.23 433.33-966.67 8.13-13.33 416.67-900.00 
(306.67+36.76) (3.9140.32) (757.99+214.24) (11.03+1.75) (619.52+179.77) 
SQW 254.09-400.00 3.45-4.31 443.65-936.67 7.64-12.67 -n 


(302.58438.66 ) (4.0240.37) — (727.674197.36) (9.69+1.23) 


ST 315.00-406.00 3.61-4.71 475.00-860.00 6.33-9.59 75.00-130.00 
(358.00424.82) (4.05+0.45)  (680.004141.39)  (8.22+1.11) (119436.39) 

IC 266.67-466.70 2.63-6.99 378.64-683.33 3.84-7.33 = 
(367.96+51.99) (4.6141.08)  (580.67487.19)  (6.16+0.99) 

LC 255.40-405.70 3.67-4.25  425.00-897.67 8.04-13.00 _ 


(304.67435.76) (3.98+0.36) (689.69+203.49) — (10.97+1.67) 


C.oligoclona 90-330 2-6 300-400 1-10 _ 
C. fracta 60-360 1-3 60-240 3-6 _ 
C. vagabunda 80-140 4-12 20-60 4-6 = 


2.3 Datasets and alignments 

The SSU and partial LSU regions of the isolated material were deposited in to the GenBank 
and the accession numbers are given in Table 1. Fifty Cladophora specimens examined, a total of 
13 ribotypes were detected (Table 1).The alignments of the two ribosomal genes, SSU and partial 
LSU, were 1 417 bp and 495 bp in length, respectively. The SSU fragment was approximately 3 
times as long as the partial LSU region but it contained about the same number of variable and 
parsimony-informative characters (SSU: 322/144; LSU: 330/155). It seems that the pairwise 
sequence divergence in the SSU was significantly lower than in the LSU (Leliaert et al, 2007). 


The alignment of SSU sequences and partial LSU sequences was 1912 bp in length, including 282 


parsimony-informative sites among 651 variable sites. The pairwise distances between studies 
haplotypes were less than 2% both in SSU and partial LSU sequences. No significant saturation 
was tested in either SSU or LSU regions, according to the /,, statistic (Xia & Xie, 2001). 
2.4 Phylogenetic analyses 

Phylogenetic analyses using BI for the combined SSU-partial LSU data (Fig.2) showed better 
performance in terms of resolution than the separate SSU and partial LSU analyses (trees not 
shown). The overall trees had similar topology for the analyzed taxa, but differed in the placement 
of a few taxa. The Cladophora haplotypes occupied a separate well-supported (0.99/0.91) position 
in most analyses sister with Cladophora vagabunda. The placement of C. sterrocladi, 
Siphonocladus tropicus and Valonia aegagropila differed among different analyses. Without 
forming a monophyletic clade, the S. tropicus and V. aegagropila clade positioned on the main 
Cladophora branch in the SSU trees (0.88) or sister with Pithophora clade in the LSU trees 


(0.82/0.78). 
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Note: A. A tree based on the SSU sequences; B. A tree based on partial LSU sequences. Other data on 


Cladophorales were downloaded from GenBank. Numbers at the nodes are BI/ML bootstrap values. Location 


abbreviations were as given in Table 1. 


Fig. 2 Hypothesized phylogenetic trees based on Bayesian analysis for Cladophora fracta and relatives 


3 Discussion 


Combining sequences from different data sets has long been questioned in phylogenetic 


analyses (Huelsenbeck et al, 1996; Leliaert et al, 2007), while the combination of multiple-gene 


data sets show better resolved and supported trees, compared with single-gene partitions (Leliaert 


et al, 2007; Boedeker et al, 2012). In the present study, the combined SSU+ LSU data were found 


to be superior to the individual partition trees, on account of the SSU fragment contained less 


variable and parsimony-informative characters and as its length was approximately 3 times as long 


as the partial LSU region. The taxonomy of genus Cladophora is problematic (Leliaert et al, 2007; 
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Ulva lactuca HQ603286 


Leliaert et al, 2009 ). Phylogenetically, all the ribotypes did not correspond with any previously 
described Cladophora fracta specimens, because we didn’t get any available sequence data of 
Cladophora fracta from genbank. However, our molecular phylogenetic analysis showed a well 
supported Cladophora clade. 

The inability to match our plants with a described taxon is confirmed by morphological data. 
Our phylogenetic analyses show that all the ribotypes formed a single clade sister with 
Cladophora vagabunda. However, C. vagabunda has distinct morphological features. For 
instances, the filament cells are cylindrical, with 80-140 um in diameter and 4-12 um in length; 
the branch lets taper to 40 um in diameter, slightly constricted at the junction with main axes; the 
apices are straight above, curved or sickle-shaped below; and the apical cells are 20-60 um in 
diameter, and 5-11 um long (Russell & Balazs, 2000), but it shares similar length/diameter ratios 
in the main axes cells with our samples. It was reported that C. vagabunda shares morphological 
similarity with C. glomerata, while differing from each other on the cell diameter and cell wall 
thickness which are easily affected by the salinity regimes. C. vagabunda and C. glomerata also 
showed an incredibly close relationship in the phylogenetic trees (Hayakawa et al, 2012), and both 
have been the most often mentioned species in eutrophic freshwaters (Whitton, 1970). C. fracta is 
separated from C. glomerata by the typically more slender filaments and the more pronounced 
tendency of the akinetes to swell into the pear-shaped structures. Considering the morphological 
characteristics, we therefore hold our species as Cladophora fracta. 

Table 2 summarizes a morphological comparison between the analyzed taxa and similar 
Cladophora species. The cell sizes were similar among/within the analyzed ribotypes; however, 
the ST samples distinguished from the XA individuals by a thinner rhizoid-like branch (75-130 
um), while XA samples exhibited the largest length/diameter ratio in the apical cells. The 
rhizoid-like branches were not observed in the SQW, JC and LC specimens. Although Cladophora 
species frequently appear in karst springs, their genetic diversity, phenological patterns and 
physiological properties have not been well investigated in these environments. In this study, 
collections from the five karst springs showed a unique morphological characteristic. Our samples 
collected from the five karst springs did not exhibit significant variation in cell dimensions, but the 
cell diameters of the individuals from the five locations were larger than those reported for the 


cells from lakes (Whitton, 1970), indicating that the habitat conditions may have an effect. A 
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similar tendency was observed in cell length and cell wall thickness. Perhaps these differences 
may have been caused by regional differences. Sinha (1968) had concluded that the number of 
nucleoli and chromocentre like bodies can be correlated with the degree of polyploidy which 
makes species larger in size. Moreover, morphological plasticity is related to diverse habitats and 
different salinity levels (Nienhuis, 1975; Hayakawa et al, 2012; Ichihara et al, 2013). Some 
frameworks, however, still lack enough sequence data to support. It is clear that additional 
sampling of Cladophora species from different salinity levels as well as karst springs will be 


needed to further clarify the diversity and plasticity of the species. 
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